Single particle Lagrangian dispersion in Bolgiano turbulence 
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Single-particle dispersion in the Bolgiano-Obukhov regime of two-dimensional turbulent convec- 
tion is investigated. Unlike dispersion in a flow displaying the classical K41 phenomenology, here, 
the leading contribution to the Lagrangian velocity fluctuations is given by the largest eddies. This 
implies a linear behavior in time for a typical velocity fluctuation in the time interval t. The contri- 
bution to the Lagrangian velocity fluctuations of local eddies (i.e. with a characteristic time of order 
i), whose space/time scalings are ruled by the Bolgiano-Obukhov theory, is thus not detectable by 
standard Lagrangian statistical observables. To disentangle contributions arising from the large 
eddies from those of local eddies, a strategy based on exit-time statistics has successfully been ex- 
ploited. Lagrangian velocity increments in Bolgiano convection thus provide a physically relevant 
example of signal with more than smooth fluctuations. 
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Understanding the statistical properties of particle 
tracers advected by turbulent flows is a fundamental 
problem in turbulent research and a key ingredient for 
the development of stochastic models of Lagrangian dis- 
persion plMEl. In the recent years there has been great 
improvements in theoretical Q, experimental |E IE 
and numerical d, understanding of Lagrangian turbu- 
lence. Most of the studies have been concentrated on the 
so-called single particle dispersion in which the statistical 
objects are the Lagrangian velocity differences following 
a single trajectory 6v(t) — v(t) — v(0). For this quan- 
tity, dimensional analysis in statistical stationary, homo- 
geneous and isotropic turbulence, predicts (Sv 2 ) ~ et 
(where e is the mean energy dissipation). The remark- 
able coincidence with diffusive behavior is at the basis of 
stochastic models of turbulent dispersion. 

In this Letter we investigate on the basis of Direct Nu- 
merical Simulations (DNS) the statistics of single parti- 
cle dispersion in two-dimensional Boussinesq convective 
turbulence. This system is characterized by an inverse 
cascade with scaling exponent in agreement with the 
Bolgiano-Obukhov theory of turbulent convection [To| . 
As a consequence, we expect a non-diffusive behavior for 
Lagrangian velocity variance. We show that a careful 
statistical analysis, based on exit time statistics, is nec- 
essary in order to observe the expected scaling. 

The two-dimensional Boussinesq turbulent convection 
is described by the following set of partial differential 
equations [TT| : 



maintaining a mean temperature profile (T(r, t)) = G-r, 
with a constant gradient G pointing in the direction of 
gravity. 



d t w + v ■ Vw 
d t T + v ■ VT 
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kAT (1) 



where T is the temperature field and co = V x v is the 
scalar vorticity, g is the gravitational acceleration, /3 is 
the thermal expansion coefficient and k and v are molec- 
ular diffusivity and viscosity. Energy in JJ| is injected by 
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FIG. I: Energy spectrum from a Direct Numerical Simula- 
tions of Boussinesq equations at resolution N = I024 2 in 
stationary conditions. The line represents the dimensional 
prediction E(k) ~ fc~ 11//5 . The inset shows the spectrum 
compensated with the prediction. 

In equation 0]) the temperature field affects the vor- 
ticity through the buoyancy force, thus providing an ex- 
ample of active tracers. At large enough values of /3, 
the buoyancy force can equilibrate the inertial terms in 
the velocity dynamics, thus providing the mechanism at 
the basis of the Bolgiano-Obukhov theory of turbulent 
convection. 

Let us briefly recall the phenomenology of two- 
dimensional turbulent convection. The balance of buoy- 
ancy and inertial terms in equation introduces the 
Bolgiano length scale Ib ^3- At small scales, r <C Ib, the 
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inertial term is larger than buoyancy force and the tem- 
perature is basically a passive scalar. At scales r ^> l B , 
buoyancy dominates and thus affects the velocity field 
in the inverse energy cascade. At variance with Navier- 
Stokes turbulence, kinetic energy is here injected at all 
scales in the inverse cascade. The energy flux at scale r 
is obtained by a dimensional argument as: 

e r = Pg ■ (v(x +r,t)T(x,t)) ~ r 4/5 (2) 

which gives the dimensional prediction for the velocity 
structure functions: 

S p (r) = ((5 r v) p ) ~ (e r r) p/3 ~ & (3) 

with ( p = 3p/5. The exponent of the second-order struc- 
ture function (£ 2 = 6/5) gives the scaling exponent for 
the energy spectrum E{k) ~ fc~ n / 5 which is shown in 
Fig. ^ It is worth remarking that detailed numerical 
simulations have shown that velocity fluctuations display 
self-similar statistics without intermittency corrections 

m 

The dimensional prediction for the statistics of La- 
grangian velocity increments Stv = v(t) — v(0) is the 
following. Considering the velocity v as the superposi- 
tion of the contributions from eddies of different sizes, 
the variation Stv over a time t will be given by the super- 
position of the variations associated to the eddies. The 
eddies at scale r, with a characteristic time r r «f will 
be decorrelated and thus give no contribution. The ed- 
dies for which r r ~ t are the first ones to be considered. 
Assuming a scaling exponent h for the velocity the char- 
acteristic turnover time time is r r ~ Ti(r/L) 1 ~' i (tj, be- 
ing characteristic time of eddies at the large scale L) and 
thus the scale of the eddies which decorrelates in a time 
t is given by r ~ L{t/TL) l ^ x ~ h ^ ■ Their contribution to 
the Lagrangian velocity fluctuation is thus estimated as 

5 t v~v L {r/L) h ~v L (t/T L ) q (4) 

where vl represents velocity fluctuations at large scales 
and q = h/(l — h). 

On top of this, the eddies at the large scales of the order 
of L have to be considered. At these scales t -C tl and 
the contribution to velocity fluctuations is differentiable, 
i.e. 5tv ~ [dtVLjt. The typical velocity fluctuation on 
the time interval t is then given by the superposition of 
two contributions: 

StV ~ T L (d t V L )(t/T L ) + V L (t/T L )<>. (5) 

At small t/rt the dominant term will be the one with 
minimum exponent: Sv ~ mm(l,q). In the framework 
of the classical K41 theory (h = 1/3, q = 1/2) the dom- 
inant contribution in JSJ is the local one which leads to 
the diffusive behavior Stv ~ t x l 2 . In the present case of 
Bolgiano convection we have q = 3/2 and thus velocity 
increments in the inertial range are dominated by the 
infrared term cx t. Therefore, a standard analysis of ve- 
locity fluctuations, i.e. Lagrangian structure functions 



Sp(t) = {(8tv) p ), is unable to disentangle the Bolgiano- 
Obukhov scaling in the Lagrangian statistics (see Fig.[2J). 

Lagrangian velocity increments in Bolgiano convection 
provide a physically relevant example of signal with more 
than smooth fluctuations. The statistical analysis of this 
kind of signal has been recently addressed on the basis of 
an exit-time statistics |13| . In the following we will dis- 
cuss the application of this approach to the present case 
and the resulting bifractal distribution of scaling expo- 
nents. 
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FIG. 2: Second order Lagrangian structure function. The 
dominant ballistic regime (Stv ~ t) is evident. 

The exit-time statistics is based on the time increments 
T(Sv) needed for a tracer to observe a change of Sv is its 
velocity 0|- Given the set of thresholds Sv n = p n Sv(0) 
(with p > 1), one computes the exit-times T n correspond- 
ing to each threshold following the tracer trajectories. 
Now, let us assume to have a signal Sv composed by 
two contributions as in (J5J. In the limit of small Sv, 
the differentiable part oc t will dominate except when 
the derivative d t VL vanishes and the local part thus be- 
comes the leading one. For a signal with 1 < q < 2, 
its first derivative is a one-dimensional self-affine signal 
with Holder exponent £ = q — 1, which thus vanishes on 
a fractal set of dimension D = 1 — £ = 2 — q. Therefore, 
the probability to observe the component 0{t q ) is equal 
to the probability to pick a point on the fractal set of 
dimension D, i.e.: 

P(T - Sv 1 '" 1 ) - T l -° - {Sv) 1 ' 1 ^ . (6) 

By using this probability for computing the average p- 
order moments of exit-time statistics one obtains the fol- 
lowing bi-fractal prediction 

(T p (Sv)) - Sv x(p) , with X (P) = minfp, -+1--) . (7) 

q q 

According to prediction (JJJ, low-order moments (p < 1) 
of the inverse statistics only see the differentiable part of 
the signal, while high-order moments (p > 1) are domi- 
nated by the local fluctuations 0(t q ). 
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FIG. 3: The spectrum E(u>) of the Lagrangian velocity incre- 
ments Sv. The line represent the prediction w~ 4 . The inset 
shows the compensated plot E(ui) /uj~ 4 . The u>~ 4 contribu- 
tion is due to the superimposed more than smooth fluctua- 
tions, while the tail of the spectrum are dominated by the 
contributions due to the non-periodicity of the signal. 



We now turn to the numerical procedure. The tur- 
bulent velocity field is obtained by direct integration of 
convective Navier-Stokes equation with a standard, 
fully-dealiased pseudospectral method in a doubly peri- 
odic square domain of resolution 1024 2 . For numerical 
convenience, the viscous term in Q has been replaced 
in our simulations by a hyperviscous term of order eight. 
Time evolution is implemented by a standard second- 
order Runge-Kutta scheme. Equations are integrated for 
about one hundred of large-scale eddy turnover times to 
reach a stationary state. Lagrangian trajectories are then 
obtained by integrating x(i) = v(x(t),t) with the veloc- 
ity at particle position obtained by linear interpolation 
from the nearest grid points. A single run follows 64000 
particles, homogeneously distributed on the integration 
domain. Velocity variations Stv and exit times T{5v) of 
single particles are recorded for times comparable to the 
integral scale. Average is performed over about 150 runs. 

In Fig. |3 we show the variance of Lagrangian velocity 
fluctuations as a function of time. The smooth behavior 
((5tv) 2 ) ~ t 2 is clearly observable making impossible the 
observation of the Bolgiano-Obukhov scaling. We remark 
that the same result is expected for any system with large 
scale domination over the local contribution, i.e. for any 
q < 1. Therefore, for these systems (with h < 1/2) the 
direct computation of Lagrangian structure functions (jSJ) 
is unable to disentangle the turbulent components in the 
velocity field. 

A first check of Bolgiano scaling in Lagrangian statis- 
tics is obtained by computing the temporal energy spec- 
trum E(u). We remind that for an infinitely differen- 
tiable signal, the smooth contribution gives an exponen- 
tial spectrum, while the S t v ~ i 3 / 2 contribution gives a 
spectrum proportional to w -4 . Figure |31 clearly shows 
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FIG. 4: Exit times (T p {Sv)) of of order p = 1,2,3 vs Sv. 
Lines represent the predicted power-law behavior JJJ with 
exponents x(p)- 

the expected ui~ 4 behavior for an intermediate range of 
frequencies. The high frequencies are dominated by the 
effects of non-periodicity of the signal. 

We now turn to the exit-time analysis. Figure |U shows 
the first moments of exit times (T p (5v)) which display a 
clear power-law scaling in the range 2 x 10~ 3 vl < Sv < 
5 x 10~ 2 vl. We were able to compute the moments up to 
p = 3 with statistical significance. The set of scaling ex- 
ponents obtained from a best fit in this range is shown in 
Fig-El The bifractal spectrum predicted by (J7J) is clearly 
reproduced. We remark that the fact that for p > 1 ex- 
ponents follows the linear behavior x(p) = (2p+ l)/3 in- 
dicates the absence of intermittency in Lagrangian statis- 
tics. This feature is a consequence of the self-similarity 
of the inverse cascade in two-dimensional Bolgiano con- 
vection. 




FIG. 5: Scaling exponents x(p) f° r the moments of the exit 
times (T p (8v)). The lines represent the bifractal prediction 
@. The error bars on the exponents have been estimated by 
evaluating differences in x(p) changing the fitting interval. 
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In conclusions, we have analyzed single-particle disper- 
sion in the two-dimensional Bolgiano turbulent convec- 
tion. A remarkable property of this turbulent regime is 
that the single-particle dispersion turns out to be dom- 
inated by the large-scale eddies whose contribution to 
velocity fluctuations behaves linearly with time. A com- 
pletely different scenario thus emerge with respect to the 
classical single-particle dispersion in a flow which dis- 
plays the standard K41 phenomenology. The Bolgiano- 
Obukhov contribution to the scaling of Lagrangian ve- 
locity increments thus results undetectable by standard 



Lagrangian structure functions. To disentangle its effect 
from the (trivial) one played by large scales, we have ex- 
ploited a method of analysis based on exit-time statistics 
for signals having more than smooth fluctuations. The 
obtained results clearly permit to identify the contribu- 
tion of " local" eddies to the Lagrangian velocity fluctua- 
tions. 
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